function y = mccormck(x)
% MCCORMCK  function
% x0=[1,...,1]
n=length(x);
y=0;
for i=1:n-1
    y=y+(-1.5*x(i)+2.5*x(i+1)+1+(x(i)-x(i+1))^2+sin(x(i)+x(i+1)));
end


